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Abstract 

N=2 supersymmetric Yang-Mills theory with flavor hypermultiplets at finite 
temperature and in the 0IS4 are studied for finite quark number density (n&) by 
a dual supergravity background with non-trivial dilaton and axion. The quarks 
and its number density rib are introduced by embedding a probe D7 brane. We 
find a critical value of the chemical potential at the limit of rib = 0, and it 
coincides with the effective quark mass given in each theory for rib = 0. At this 
point, a transition of the D7 embedding configurations occurs between their two 
typical ones. The phase diagrams of this transition are shown in the plane of 
chemical potential versus temperature and cosmological constant for YM theory 
at finite temperature and in dS4 respectively. In this phase transition, the order 
parameter is considered as rib- This result seems to be reasonable since both 
theories are in the quark deconfinement phase. 
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1 Introduction 



In the context of the holography [H |2] , the dynamical properties of flavor quarks have 
been studies by embedding D7 brane(s) as probe. [3j HI El El [7J El [9j [101 IH] • Recently, 
the research in this direction has been extended by introducing quark number density. 
Such researches have been performed at finite temperature for D3/D7 [TBI [T4"l [T5] and 
D4/D8 [12] models, and some kinds of phase transitions are reported. 

In D3/D7 model, the high temperature gauge theory is represented by the AdS 
Schwartzchild configuration, which describes chiral symmetricB and deconfinement phase 
for zero quark number . A horizon appears in the bulk configuration in this model, 
then two types of D7 embeddings are possible, the "Minkowski embedding" , in which 
the D7 brane is off the horizon, and "black hole embedding" (BH), where the IR end 
point of the D7 brane touches on the horizon. 

For the case of non-zero n&, however, the authors of p3] have pointed out that 
the Minkowski embedding is unphysical due to the fact that the inner electric field 
remains finite at the IR end point of D7 brane, and it can not go nowhere from that 
point without any appropriate source like a baryon. Therefore, the physical embedding 
is restricted to the BH type, and they found that all the quark mass range could be 
covered by the BH embedding. This implies that the BH embedding is enough for 
finite rib case. 

However, we show here that a finite region of chemical potential (0 < ix < no) can 
not be covered by finite rib and BH embedding, and this region is described rather by 
the physical Minkowski embedding with rib = 0. And the critical point /io(= f^oiT)), 
which depends on the temperature T, is shown to be equivalent to the effective quark 
mass (rh q (= rh q (T))). Then rib appears when ii exceeds rn q . This result implies that 
rib describes the number of quarks moving freely with the mass rh q in the thermal 
YM medium. This interpretation is very reasonable since the chemical potential is 
examined here for AdS Schwartzchild background which is dual to the gauge theory in 
the high temperature deconfinement phase. Then it would be an open problem here 
how rib could be related to the baryon number density of a physical baryon. And the 
curve fio{T) in the ii-T plane indicates the critical line which separates two phases, 
Minkowski embedding with rib = and BH embedding with rib 7^ 0. In passing this 
critical curve, the chemical potential does not jump but the electric potential and 
profile function of D7 embedding show a jump near the horizon. The order parameter 
of this transition is interpreted as rib oc (^ r '^ r ), where \I/ denotes the quark field. 

Similar situation to the finite temperature gauge theory is seen also for the gauge 
theory in the de Sitter space time (dS&) [2TJ [22] . This theory is characterized by the 
4d cosmological constant (A), and the configuration of the gravity dual has a horizon 
as in the high temperature theory. Then both the "Minkowski" and BH D7 brane 
embeddings are possible. The dynamical role of A is similar to the temperature so 
far. In fact, we find for finite A that the chiral symmetry is restored and the quarks 
are deconfined. Then it is useful to study the meaning of the "baryon number" rib 
also in the dS^ gauge theory. Actually, we find a similar critical curve /xo(A) in the 

* Here and hereafter we use the chiral symmetry for the state with massless quark or a hypermul- 
tiplet, and the breaking of this symmetry means the spontaneous mass generation for the quark. 
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/x-A plane and the same kind of phase transition on this critical curve. We also show 
/x (A) = fh q (X), 

where rh q (\) represents the effective mass of a quark moving freely in dS4 space. 

Thus throughout the gauge theory at finite temperature and in dS^ we find that 
the chemical potential remains finite as /xo in the limit of rib — in a definite region 
of T or A. In other words, there is an area of n& = and finite /x(< /xo) in /x-T and 
/x-A planes. In this area, the D7 brane is embedded by Minkowski type, but the finite 
chemical potential plays no dynamical role since it disappears in the brane action since 
the electric potential is given by a constant /x. However, when it exceeds a critical 
value /Xo, ^fe(oc (\Ed\l/)) appears as a non-zero value and the D7 profile is changed to 
BH type. This phenomenon could be interpreted as a phase transition with the order 
parameter (\Ed\l/). Our purpose is to give such phase diagrams and to discuss this 
transition observed in the deconfinement phase. 

In section 2, we give the setting of the bulk background and the D7 brane action 
and the equations of motion of embedding are given in section 3. In the section 4, the 
chemical potential and its relation to the phase transition is given, and phase diagrams 
are given in the section 5 and 6. The summary is given in the final section. 



2 Background geometry 

We start from lOd IIB model retaining the dilaton <3>, axion x an d selfdual five form 
field strength F^y Under the Freund-Rubin ansatz for F( 5 y F^...^ = —\fk/2 e M1 ... M5 
[T6j [TTj . and for the lOd metric as M 5 x S 5 or ds 2 = gMN<lx M dx N + g^dx l dx^ we find 
the solution. The five dimensional M 5 part of the solution is obtained by solving the 
following reduced 5d action, 

S = ^ / d 5 xy/=g- (R+3A- \m 2 + \e 2 *{d X ) 2 ) , (1) 
which is written in the string frame and taking a' — g s — 1. 

Two deconfinement solutions 

The following solutions are obtained under the ansatz, 

X = -e~* + Xo , (2) 

which is necessary to obtain a supersymmetric solutions. However, the supersymmetry 
is broken in the following two solutions due to the cosmological constant A and the 
temperature T. 

(2) High temperature phase: First, we consider the high temperature gauge 
theory given in [20J, and it is written as 

ds^Q = G MPfdX^ dX^ 

= e* /2 {-f{r)de + (dxr) + J^T^dr 2 + , (3) 
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i( i * ■' i 

/ l-(r T /r) 4 ' 



e*= 1 + 4log( t ^J ) , (4) 



/(O = yjl - (y) 4 , (5) 

where where M, N = ~ 9 and i? = = (47T./V) 1 / 4 . The temperature T is 

related to the parameter tt as ry = ttR 2 T, and g represent the VEV of gauge fields 
condensate [10], q oc (F 2 ). 

(2) dS 4 : The second solution is for the dS 4 solution 



dsl = e*/ 2 (^A 2 (-rft 2 + a{t) 2 {dx*) 2 ) + ^dr 2 + R'dQ 2 ] , (6) 



. q 1 - (r /r) 2 /3 

' r 4 (1 - (r /r) 2 ) 3 ' 1 J 

A = l-(^) 2 , a(t)=e 2 ^ (8) 
r 

The horizon is denoted by r , which is related to the 4d cosmo logical constant A as 

A = 4 (9) 

In this configuration, the four dimensional boundary represents the jV=4 SYM theory 
in the de Sitter space or in the inflational universe characterized by the 4d cosmological 
constant A. 



3 D7 brane action and equation of motion: 

The embedding of the brane is performed by solving the equations of motion for the 
fields on the brane. The brane action for the D7-probe is given as 

Sdy = -T 7 J d s £ (Vy-det (G ab + 2ira'F ab ) - ^ ;: "\ 1,,..,^ 

+ { ^p^r 7 Jp[C^}AFAF, (10) 

where F ab = d a A b -d b A a . Q ah = d^X M d^X N G MN (a, b = ~ 7, and M, = ~ 9) 
and r 7 = [(27r) 7 (? s a' 4 ] _1 represent the induced metric and the tension of D7 brane 
respectively. And PfC^ 4 -*] denotes the pullback of a bulk four form potential, 

C (4) = (j^dx A dx 1 A dx 2 A dx 3 ^j . (11) 

The eight form potential Ai v ..i 8 , which is the Hodge dual to the axion, couples to the 
D7 brane minimally. In terms of the Hodge dual field strength, F( 9 ) = cL4( 8 ) [T8] . 
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the potential A(s) is obtained. In this brane action, there are two scalar fields, which 
are denoted here as X 8 and X 9 , corresponding to the two coordinates transverse to 
the brane. Since the brane action is rotational symmetric in the plane of these two 
coordinate, the embedding is given by the solution of one of the scalar fields, and it is 
denoted by w(= X 8 ) given below. 

In addition, in the present case, we must solve for the Wick rotated time-component 
(Aq) of the U(l) vector fields on the brane in order to introduce the chemical potential 
(fi) of the quark. 

For finite temperature : 

In the case of finite temperature, the D7 brane is embedded as follows [20]. The 
six dimensional part of the above metric is rewritten as, 

' R2 -dr 2 + R 2 dn 2 5 = ^ (dp 2 + p 2 dn 2 3 + (dX s ) 2 + (dX 9 ) 2 ) , (12) 



f 2 (r) r 2 5 U 2 



U(r) = exp I / , t = Jl±M. (13) 



rJl — (rx/rY 
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Here U is normalized as U = r for ry = 0, and U 2 = p 2 + (X 8 ) 2 + (X 9 ) 2 . Then we 
obtain the induced metric for D7 brane, 



ds 2 8 = e*' 2 \ ^ (-f{r)dt 2 + (dxy) 



, R2 

^ ((i + (oyr + (ay) 2 )d P 2 + P 2 dnt) } , (i4) 

where w 8 (p) and w 9 (p) are the scalars which determine the position of D7 brane. They 
are solved under the ansatz that they depend on only p. Further we can set w 9 = 
and w 8 = w(p) without loss of generality due to the rotational invariance in X s -X 9 
plane. By substituting the above metric in the DBI action and retaining A (p), we 
obtain 



where F pt = 27ia'F pt , F pt = d p A and the last term in the parenthesis denotes the eight 
form part. 

The equations of motion are solved as follows. First, we solve the equation of motion 
for Aq by introducing the quark number density rib as an integral constant as follows, 

M*?*** -n b (16) 
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From this, we obtain the following form of U(l) electric field, 



U \n 2 + ffi(r/U)*e* 



(17) 



Here we see the property of F p t t T near the horizon r = tt- At this limit, we have 



Fpt,T ~ n b (l ) 



;is) 



for both q = and q 7^ 0. Then F pt ^ = at the horizon, while it is finite at p = 
for rib 7^ 0. So the only the BH embedding is physical for rib 7^ 0, and the Minkowski 
embeddings are considered only for n b = 0. 

Finally, then, the equation of motion for w is given by 

2(1-/) 



w" w' w — pit;' 



1 + w' 2 ' p £/ 2 / 



2 _ ^ 
r 4 



— Gj' 

where Gt is defined as 



U 2 f 

3w' (w'p — w) 

~ P Tpf~ 



s(i-f)f + % 



G 7 



TT nl + fP{r/U)*e* ' 
Of course, this equation reduces to the one given in [20] when rib — ► 0. 



(19) 



(20) 



For dS A : 

In this case, the extra six dimensional part of the above metric ([6]) is rewritten as, 



d2 d2 

^-dr 2 + R 2 dQ 2 5 = tL 



(dp 2 + P 2 dn 2 3 + (dx 8 ) 2 + (dx 9 ) 2 



(21) 



where r = p + (X ) + (X ) . And we obtain the induced metric for D7 brane, 



dsl = e^ 2 I T —A 2 (-dt 2 + a(t) 2 (dx 1 ) 2 ) + 



R 2 



((1 + (d p w) 2 )d P 2 + P 2 dn 



(22) 



where we set as X 9 = and X 8 = w(p) as mentioned above. Then retaining the profile 
function w(p) and Ao(p), we arrive at the following D7 brane action, 



2 — - c» 



A 2 e* 



(23) 



where F, 



pt 



2-Kot'Fpt, Fpt = dpAo and the eight form part is given as 
C 8 (r) = fdr' A\r')d r , (exp(<f>(r'))) = 



(24) 



The equations of motion are solved as above. First, we solve the equation of motion 
for A by introducing the quark number density n& as an integral constant, 



p 3 A 2 F pt 



l + (w 



l\2 



F 2 
A 2 e* 



n b 



(25) 



From this, we obtain the following form of U(l) electric field, 

F pt = n b Ae*/\ 



w 



A 2 



\n% + p 6 A 6 e 



(26) 



Then we get the equation of w by substituting this into the one obtained from D7 
action by the variational principle. The final form is written as, 



w 



VI + W 2 1 1 + w' 2 



+ w' 



- + ($ + 4 log A)' - G - + (3 log A + - 



+ 




w 



p + ww' 



w 



I 2 



;$ + 41ogA)'-G(| + 31ogA)' 



where prime denotes the derivative with respect to p, and G is defined as 



G 



nt, 



(27) 



(28) 



When we take the limit n b — > at non-zero p, the above equation fT27|) reduces to the 
one given in |22j . By giving n b as a parameter, we can solve this equation and find the 
profile of the embedded D7 brane. Then we find F pt simultaneously through 



4 Chemical potential and phase transition 

From the equations given above, we firstly solve for w with the boundary condition, 
w(oo) = 2'Kot'm q) where m q denotes the current quark mass and we set as 27ra' = 1 
hereafter. Then we can read from the asymptotic form of this solution the chiral 
condensate (ty^, where \l/ denotes the quark field. Then we solve for the gauge 
potential A (p) by using this solved w. 

In this case, we can assume its asymptotic form at large p as in the case of w. In 
the context of AdS/CFT, it would be written in terms of the chemical potential p and 
the quark (or baryon) number density as, 

A)(p) = ^-^ + "--, (29) 
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We notice 

n b oc = (#7°*) . (30) 

This would represent the baryon number density when the theory is in the confinement 
phase, but we call this as quark number here since the present model describes the 
deconfinement phase. 

In the present case, we firstly give n&, then solve for w(p, n h ). And by using this 
solution w, we obtain the chemical potential p = A (oo) from the following formula 

roo 

p = / dpF p t = A (oo) - A (p rn ) = p - A (p m ) , (31) 

J Pm 

where p m is the minimum value of p and we remind F pt = 2ira'F pt . 

In order to make clear our viewpoint, we firstly restrict the analysis to the case 
of large m q where w(p, = 0) is the Minkowski embedding solution. For n& = 0, 
F pt = 0, then the Minkowski embedding is physical. This does not mean p = 0, but 
it implies A = p where p is conatant. In this case, p = and A (p m ) = A (0) = p. 
So we can consider two types of solutions for A depending on the solution of w, (A) 
Minkowski type with = and (B) black hole type with ^ 0, in the case of large 
m q as given below. 

For the BH embedding solutions with n& 7^ 0, the infrared boundary value A (p m ) 
is estimated by using Eqs. (TlTl) and (1261) . For enough small rib, pm, is very small and 
the solution w very rapidly increases from the point p = p m . It soon arrives at the 
value near the maximum value w(oo) at about p = p m + e, e <C 1. In other words, w 
is approximated by a step function then w' is done by the 5 function. As a result for 
1 > n b > p z m , we obtain 



/ dp ~ / dp Ae^w' = / dr Ae* /2 , (32) 

J Pm J Pm Jro 

for dS^ with q = 0, and 

/ dp Fp t , T ~ / dp f-e^ 2 w' = / dr e*/ 2 , (33) 

J Pm. J Pm. U J TT 



for finite temperature for any q. Here we notice the dilaton denoted by $ in the above 
two equations are different from each other. Their explicit forms are given in the section 
two. For the case of 0LS4, however, the BH embedding is obtained only for q = since 
Fpt diverges at p m for q ^ 0. Meanwhile such a situation can not be seen for the finite 
temperature phase, so we can consider any value of q in this case. On these points, we 
discuss again in the below. 

In the limit of p m = 0, the most right hand side of (I3"2"|) and (I3U1) are approximated 

as 

/ drAe^ 2 = m q , / dr e^ 2 = m q: T . (34) 

They are equivalent to the effective quark mass given in [22] and [20] defined from the 
Wilson line in dS± and at finite temperature respectively for the Minkowski embedding 



7 



case with F pt = 00 In both cases, the quarks are deconfined and moves in the cor- 
responding gauge field medium with this effective mass. On the other hand, the left 
hand side of (1321) and (|33|) , by its definition, is equivalent to p — A (p m ) at p m — > 0. 
If we consider the chemical potential should be the effective quark mass in the limit 
rib 0, then we should take as 

A (p m )=0, (35) 

and we find in the limit of rib — > 

/2(n& — > 0) = /i(n& — > 0) = fiQ = rh q or m (?j t- (36) 

We must notice here that this analysis is performed for nj, ^ even if the limit of 
rib —> is taken, and the embedding profile is assumed to be the BH type. 

The important point to be noticed here is that the value of po is finite. Then there 
is a region, < p < po and rib = 0, for finite A or T in both cases. In these regions, 
the Minkowski embedding is allowed, then we can consider as mentioned above the 
following phase 

Phase (A) n b = : {p m = , A (0) = //} 
A = fi = const. . 

On the other hand for fiQ < /i, the following BH embedding phase (Phase (B)) is 
realized, 

Phase (B) n b > : {p m > , A (p m ) = 0} 

... rib , 
Ao(p) = p - ^2 H at lar S e P ■ 

And the critical curves separating (A) and (B) are given by rh q (X) and fh q> t{T) for 
each model. 

For this definition of two phases, we show the phase diagrams in the next section 
including numerical analyses. We should notice here that, in this phase transition, from 
pahse (B) to (A), A jumps from A (p m ^ 0) = to A (0) = p at the limit of p m = 0. 
Simultaneously, the embedding profile w changes from the BH to the Minkowski type. 
Actually from the left equations of the above approximate formula (|32|) and (|33|) . we 
find A' Q ~ w', then we can understand the jump of A corresponds to the jump of w 
near the horizon. 

5 Phase transition at high temperature phase 

Here we firstly solve the embedding equations by introducing the chemical potential 
and noticing that any Minkowski embedding solution is unphysical for rib ^ 0. As 
mentioned above, the Minkowski embedding is considered only for nj = 0. In this case, 

t The authors of [14j discuss on this point related to the string energy from a similar viewpoint. 
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F pt = at any point, then we find physical Minkowski embeddings but there is no 
dynamical role of chemical potential since it vanishes completely from the action. 

As stated above, two types of embeddings, (A) Minkowski and (B) BH embeddings, 
are possible, and they are discriminated by the form of A Q and w. And the critical 
curve separating the two regions is given by po = fh q . Before giving this critical curve, 
we discuss another transition which has been observed in 



(i) Phase transition within (B) at small p: 

In the high temperature model, the topology changing phase transition is observed 
for q = and n& = 0. This transition is characterized by the jump of w(p m ), where p m 
represents the minimum point of p. In this sense, this transition is the first order since 
the D7 energy jumps at this point due to the different configuration of w. This kind 
of transition is also observed for finite rib between the same BH embedding. So this is 
not the topology changing phase transition in the case of rib ^ 0. 

At the critical temperature T\, we find three different solutions at very small rib 
for the same quark mass, m q = w(oo). They are shown in the Fig. [TJ This is firstly 
observed in [13]. In [13], w is assigned as the radial coordinate in X 8 -X 9 plane, but 
this is not important. We set as w = X 8 , and we obtain almost the same result with 

The solutions are separated to two categories whether the limit value p\ nb=0 (= p ) 
is finite (sol. (a)) or zero (sol. (b)). We notice here that the two solutions for the case 
of solution (b) are overlapped in the figure due to too small n b . Then, this gives a gap 
of po(T) at this temperature 7\ as seen in the phase diagram in p-T plane given in the 
Fig. [31 The solution (a) provides finite po, and we observe 

for this solution. This fact implies that this solution is unstable as pointed out in |14j . 
And this is seen for the very small p region of < p < p%, where pi depends on m q and 
is very small. As shown below, in the dS^ model, there is no such a transition point at 
finite rib so the gap is not seen. These points are assured in the phase diagrams given 
below. 

Finally, we comment on the relation between the transition discussed here in terms 
of the Fig{T] and the one observed in the case of rib — 0. The solution (a) in the 
FigJT] transits to the Minkowski solution and the solution (b) remains as BH solution 
in the limit of rib — > 0. This situation of the two BH configurations is schematically 
depicted in the FigJH In this sense, the transition considered here can be regarded as 
the limitting case of the topology changing phase transition which is seen at rib — 0. 

(ii) Phase transition from (B) to (A): 

We notice that Po{T) is zero for T > 7\ and it becomes finite in the low temperature 
side T < T\. (See Fig. [3J ) The situation is here considered from the viewpoint of grand 
canonical ensemble, where the chemical potential p is regarded as a parameter of the 
theory. In this sense, in the region of < p < po, the solutions (A) (Minkowski 
embedding) with rib = are found. When p exceeds po, p becomes active and could 
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w(p) w(p) 




Fig. 2: Embedding solutions near the transition point for q = 0, m q = 1.30916. The 
left are for T = 0.4, and n b =0, 0.000187, 0.00184, 0.0181, 0.217 from the above. The 
right is for T = 0.8, and n b =0, 5, 10. 
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generate a finite rib((x (w^)) to form the solution (B) (BH embedding). Then, for 
//o < £t, the brane profile is changed to the BH embedding solution when its D7 energy- 
density becomes smaller than the one of the Minkowski embedding with the same p but 
with rib = 0. We show in the below that this phenomenon really occurs at /i = po(T). 

This implies a topology change of the D7 embedded configuration in the bulk 10 
dimension, and, on the other hand, it is observed as the generation of (^Vl/) in the 4d 
gauge theory. The critical point po(T) is determined for fixed m q = w(oo) and T by 
taking the limit of rib = from the side of rib 7^ (in the phase (B)). 

The typical solutions of w(p) near this transition point for (T,m q ) = (0.4, 1.3) are 
shown for various p in the left of the Fig. [21 where one Minkowski solution with rib = 
and four BH solutions with finite values rib are shown. For the latter BH solutions, we 
obtain the corresponding values of p(> po), where po at this temperature (T = 0.4) is 
read from the Fig. [3] as about p Q = 0.14. We notice that we have also the Minkowski 
solution for those value of p with rib = 0. Therefore, we must compare the D7 energy 
of each BH embedding solution and the corresponding Minkowski solution at the same 
value of p in order to see which solution has lower energy. The D7 energy E^ 7 is defined 
as 

SEr = -T7 J <?£ e 3 1/2 E T , (37) 
then for the Minkowski embedding of q = we obtain 

El = f dp p 3 (^Vyl + K) 2 , (38) 

and for BH solution 

E T BH = j2 dp p 3 {^f fe^{l + {w'Y){l-G T ) , (39) 

where Gt is given in ( |20l) . Each energy diverges, so we calculate the difference of E\. { 
and E^ H) 

AE T = E T M - E T BH 

to cancel the divergence coming from large p integration. At large p, the solutions of 
Mincowski and BH embeddings approaches to the same function, then this method of 
the regularization works well. And we could see Ej^ > E BH for any case, then the 
transition occurs at po between the Minkowski solution and the one of the BH with 
the smallest rib- In other words, the curve Po{T) obtained in this way gives the critical 
curve which separates two phases (A) and (B), which are defined above. 

In order to assure the statement, po(T) = for T > Ti, given above, we also studied 
at T = 0.8 where all the solutions are BH type even if rib = 0. The solutions are shown 
in the right hand side of Fig. El In this case, we could see that Eq H decreases with 
increasing rib. So no Minkowski embedding appears as a stable state. 



The resultant curve of po{T) is shown in the Fig. [3l Near T = T\ ~ 0.45, a 
small flat region is seen. This point is the transition point shown in (i) above. The 
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value of /i jumps from to a small but finite fit. Then, for ii\ < ji, Ho(T) increases 
monotonically with decreasing T, and arrives at m q = w(oo) in the limit of T = 
0. However, this diagram has its meaning only for T > T c , where T c denotes the 
deconfinement /confinement transition temperature since the present model describes 
only the deconfinement phase except for T = and q > 0. While we do not discuss on 
this phase transition here, 
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.2 
. 1 
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Fig. 3: Phase diagram in T-/x plane for R=l,m g = 1.30916. The dots represent the 
effective quark mass rh q given by the last equation (1341 . 



it is shown up to T = 0, from theoretical viewpoint. 

The values of the effective quark mass rh q given by the last equations of (I3"4~|) are 
shown in the Fig. [3] by several dots, and we can see they precisely coincide with fio(T) 
which are obtained as above. 

At low temperature T < T c , we should consider a model in the confinement phase. 
In our model, it is realized for T = and q > 0. In this case, we know that fh q diverges 
[TU] and this means that we can not observe an isolated quark. In other words, the 
theory is in the confinement phase. But we always find a finite rh q for T > even if 
q > in the present model, so our model describes the deconfinement phase so far as 
T is finite. Then /io(0) is infinite for q > 0, but the behavior of ^o{T) for T > is 
qualitatively similar to the one of q = 0. Then, the qualitative property is well described 
by fio(T) of q = 0. An important point is that the chemical potential is inactive for 
fi < fio, and there is a threshold where it becomes dynamical and to generate non-zero 
rib oc (^f^\. This result implies that the chemical potential becomes active when it 
exceeds the effective quark mass since the gauge theory is in the quark deconfinement 
phase and there is no bound state of quarks. In order to make clear this point we 
consider the relation with a possible bound state. 

As another possibility, we can consider /i as the threshold of some bound state of 
quarks. In the present case, it would be the baryon with the mass Mb, then we will 
have fi = M B /N C [TUJ . In order to prove this conjecture, we must include the baryon 
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Fig. 4: The extended small \x region of the Fig. [3] is shown. At first, the two BH 
embeddings appear on the horizontal line at (T, /x, n b ) = (0.450158, 0.0223708, 0.00699) 
in a degenerated form. Then for decreasing n&, they run on the line in the opposite di- 
rection shown by the arrows and approache to the point (a) and (b), which corresponds 
to the two BH configurations shown in the Fig{T] 



in the model by for example introducing the D5 brane. In this paper, this point is an 
open problem and we will examine it in the near future. 

In [19], the critical value of the isospin chemical potential is shown to be equal to 
the half of the meson mass in the case of quenched approximation. So this point could 
be examined by extending the model to non-Abelian DBI action [23]. In a related 
direction, an interesting analysis in the Higgs branch has been seen [24J. 



6 Phase transition in dS± 

In this case, the bulk configuration has a horizon at r = tq which determines the 4d 
cosmological constant A as A = 4rg/i? 4 . Then the both the Minkowski BH embeddings 
are seen, and we can expect the change of profiles as in the finite temperature case. 
Actually, a kind of phase transition is seen in the case without the chemical potential 
but with gauge field condensate [22] . Here we are solving the embedding equations by 
introducing the chemical potential. 

In this case, we should notice the following points, (i) For the Minkowski embed- 
ding, F pt should vanish at the end point p = 0. (ii) As for the BH embeddings, the end 
point of the brane is at the horizon r = r Q . Near this point, we find e* oc q(l — r /r)~ 3 
then 

F pt ocg^Vl-ro/r)" 1 / 2 
which implies F pt = oo at r = tq. In this case, we can not find any black hole 
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embeddings. As a result, the black hole embedding solution can not be obtained for 
finite q (finite gauge condensate (F^ u )) and rib(^ 0). 
Meanwhile, for q = and n& 7^ 0, 

e* = 1 , F pi oc (1 - r /r) -> 

in the limit r = ?~o, and this is independent of rib- Then, in the case of q — 0, we find 
-F pt = on the horizon. It shows that the black hole embedding is physical in this case, 
and we actually could find such solutions. 

Therefore, we restrict hereafter to the case of q = for dS^ model in order to allow 
the physical embeddings of n b 7^ the black hole embedding, since it is the 

only possible embedding in the case of n& 7^ 0. On the other hand, as mentioned above, 
the Minkowski embedding is allowed only for rib — (phase (A)). But, we remind that 
this does not necessarily mean Aq = 0. For this embedding, we can consider a finite 
chemical potential, Aq — // where \i is finite but rib = 0. In this case, F pt = at any 
point, then we find physical Minkowski embeddings but there is no dynamical role of 
chemical potential since it vanishes completely from the action. As in the above case, 
the phase (A) is realized for < \i < Po(X), and /io(A) represent the critical curve of 
the phase transition from phase (A) to (B). And this curve is obtained according to 
the method given for the finite temperature case. 

Fixing m q and A, the typical solutions of w(p) for q = near this transition point 
are shown in the Fig. [5j 



w(p) w(p) 

3.5r 2.5r 




Fig. 5: Embedding solutions near the transition point for q = 0, m q = 2.94966. The 
left are for A = 4, and n b =0, 0.000124, 0.00124, 0.0133 and 0.0971 from the above. 
The right is for A = 6, and n b =0, 0.001, 0.01. 



In the left of the Fig. \5[ four BH solutions have finite values of rib as shown in the 
figure caption. They have also the corresponding values of // > /i , where /i is shown 
in the Fig. [61 For those values of n, we have also the Minkowski solution given in the 
Fig. [5] with A = p. Then, we compare the D7 energy of each BH embedding solution 
and the corresponding Minkowski solution at the same value of \x in order to see which 
solution has lower energy. The D7 energy for q = which is defined as [22] 

Sd7 = -t 7 Jd 7 ^a(t) 3 E dS4 (40) 
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then for Minkowski embedding solution, 

r°° i 

E^ = I dp pWv/l + K) 2 , 

and for BH solution, 

POO > 

= / dp pWy(l + K) 2 )(l-G) , 



(41) 



(42) 



where G is given by fl28l) . But as in the case of finite temperature, _Eds 4 is also divergent. 
Although a way to reguralization is given in [22] , we estimate the difference, 



AE = - E™ 

as in the previous section at the same quark mass m q and A since AE should be finite. 



A 



6 








5 








4 








3 




n b ^0 




2 


n b = 






1 










0.5 1 


1.5 2 


2 . 5 m q 3 


-1 









Fig. 6: The value of /io for mq=2. 94966 and R=l. The dots represent the effective 
quark mass calculated by (1541) for m q = 2.94966. 



According to this procedure, we find that all the above black hole solutions give 
lower energy than the one of the Minkowski embedding, and the energy decreases with 
increasing n^. Then the critical value po at this A is obtained as the limiting value of p 
obtained by approaching to n& = from BH embedding side. For A > 5, this limiting 
value is obtained as zero. Then the all the embeddings are the BH type even if rib — 
as shown by the right one of the Fig. [5j Performing this procedure for other points of 
A < 5, we find the critical curve in the X-p plane. The results are shown in the Fig. [6j 
As in the high temperature case, we can see the equivalence of p and the effective 
quark mass fh q obtained from the formula (|32l) . Then this result is interpreted as the 
common property of the deconfinement phase gauge theory with quarks 
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7 Summary 



Here the role of the chemical potential of the quark is studied in terms of the holographic 
gauge theory at finite temperature and in dS^. Quarks and the chemical potential are 
introduced by embedding the probe D7 branes in the bulk configurations corresponding 
to the considering gauge theories. In both bulk configurations, there is a horizon. Then 
two embedding forms, Minkowski and black hole, are possible. The chemical potential 
p and the number density n b are introduced through the Wick rotated time component 
of U(l) vector field, A , on the D7 brane. 

In order to embed the D7 brane, the profile function w(p) and A (p) are solved by 
the variational principle of the D7 action. The equations of motion are solved firstly 
for w by giving n b , then we obtain p. This method correspond to study the system 
in the micro-canonical ensemble. In this approach, we find a finite p(= po) in the 
limit of n b = 0. This implies the existence of a parameter region of < p < p and 
n b = 0, where the D7 brane can not be embedded by the BH embedding. So we need 
Minkowski embedding in this region in order to introduce quarks. 

Actually, in this region the Minkowski embedding is physical since n b = and we 
obtain the electric potential Aq = p, then Aq is a constant with respect to p. Thus, 
in this case, we can study the system from the viewpoint of grand canonical ensemble, 
where p is the given parameter. So, by varying the chemical potential from zero to 
large p, we find that p is inactive below p and (^^) appears for p < p and increases 
with p. 

At the same time with the appearance of the profile of the D7 embedding is 

changed to the BH form. This is assured by comparing the D7 energies of Minkowski 
and BH embeddings with the same p and m q at fixed T (or A). The energy of the BH 
embedding with finite n b is always smaller than the one of Minkowski embedding. This 
is therefore a kind of phase transition with the order parameter In the bulk, 

this is seen as the topology changing of the D7 embedded configuration. We assured 
this phase transition in the two different theories which are in the quark deconfinement 
phase. In both cases, the critical point po is identified with the effective quark mass 
fh q . Thus, this result is consistent with the picture that the quark in the deconfinement 
phase moves freely with this effective mass. 

Then the transition stated above is common to the two deconfinement model. But 
we observe some differences between the two models. In both cases, we introduce a 
parameter q which is representing the gauge condensate (F^ v ). And this parameter is 
responsible to the quark confinement. Actually, for A = and T = 0, the two theory 
are equivalent and the theory is in the confinement phase. On the other hand, for 
T > or A > 0, both theories changes to the deconfinement phase. However, when n b 
is added, while the BH embedding is possible in the finite temperature theory for any 
value of q > 0, it is forbidden for dS 4 . 

Another difference is a phase transition which is seen only for the finite temperature 
case near very small p, where a embedding configuration change is seen between the 
same BH type configurations. This is not seen in the case of dS 4 . 



16 



The interesting case would be the model in the confining phase, where the effective 
quark mass diverges, then /xq diverges. In other words, we can not find finite 
at any /j. The only case, where (^^) is seen in the confinement phase would be in a 
model in which the baryon is introduced. In this case, we will find the transition point 
at /Jo = Mb/N c through the same analysis, where Mb is the lowest baryon mass. On 
this point we will discuss in the future. 
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